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ABSTRACT 

Precession in an accretion-powered pulsar is expected to produce characteristic vari- 
ations in the pulse properties. Assuming surface intensity maps with one and two 
hotspots, we compute theoretically the periodic modulation of the mean flux, pulse- 
phase residuals and fractional amplitudes of the first and second harmonic of the pulse 
profiles. These quantities are characterised in terms of their relative precession phase 
offsets. We then search for these signatures in 37 days of X-ray timing data from the ac- 
creting millisecond pulsar XTE J1814— 338. We analyse a 12.2-d modulation observed 
previously and show that it is consistent with a freely precessing neutron star only if 
the inclination angle is < 0.1°, an a priori unlikely orientation. We conclude that if 
the observed flux variations are due to precession, our model incompletely describes 
the relative precession phase offsets (e.g. the surface intensity map is over-simplified). 
We are still able to place an upper limit on e of 3.0 x 1CP 9 independently of our model, 
and estimate the phase- independent tilt angle 6 to lie roughly between 5° and 10°. On 
the other hand, if the observed flux variations are not due to precession, the detected 
signal serves as a firm upper limit for any underlying precession signal. We then place 
an upper limit on the product ecos9 of 9.9 x 10~ 10 . The first scenario translates 
into a maximum gravitational wave strain of 10 -27 from XTE J1814— 338 (assuming a 
distance of 8 kpc), and a corresponding signal-to-noise ratio of ^ 10~ 3 (for a 120 day 
integration time) for the advanced LIGO ground-based gravitational wave detector. 

Key words: gravitational waves — pulsars: general — pulsars: individual (XTE 
J1814-338) — stars: neutron — stars: rotation — X-ray: binaries 



1 INTRODUCTION 

Accreting millisecond pulsars (AMSPs) are a subset of neu- 
tron stars in low-mass X-ray binaries (LMXBs) that exhibit 
persistent X-ray pulsations with periods below 10 ms. In 
the standard recycling scenario, AMSPs are the evolution- 
ary link bet ween LMXBs and nonaccreting, radio millisec- 
ond p ulsars (lAlpar et al.l ll982; Radhak rishnan fc SrinivasarJ 
Eight AMSP s have been discovered at the time of 



writing dWiinands 


l2004l;lMorean et al.ll2005l:lGallowavll2007l: 


IKrimm et al.l 2007 


)• 



Most AMSPs are X-ray transients. Once every few 
years, they emerge from quiescence and become detectable 
during an outburst lasting several we eks. Th e outb urst is 
attributed to enhanced accretion [e.g. lLasotal (|2001T )]. fun- 
nelled onto a small number of hotspots on the star. Lit- 
tle is know n about the shap e , posi ti on, or number of these 
hotsp ots |Romanova et al. | 12004 iKulkarni fc Romanoval 
l2005h . but they do give rise to detectable X-ray pulsations, 
from which the spin period and orbital parameters can be de- 



termined. During an outburst, surface thermonuclear burn- 
ing also causes type I X-ray bursts, which last a few min- 
utes and occur on average once every few days. Type I 
X-ray bursts have been observed in t hree AMSPs to date: 
SAX J1808.4-3658, X TE J1814-338 JWiinandsll2006l '). and 
HETE J1900. 1-2455 (|Vanderspek et al.ll2005l )~In the burst 
tails, a small component of the X-ray flux (~ 15% for XTE 
J1814— 338) oscillates at the spin frequency. 

AMSPs are expecte d to be rela t ively powerful grav- 
itational wave sources (I Watts et all 12008). The fastest, 
ICR J00291+5934 (jEckert et al.l |2004 ). spins at Q.*/2n = 
599 Hz, well below the theoretical breakup fr equency for 
most nuclear equat ions of state (~ 1.5 kHz) (|Cook et al.l 
1 1 994 iBildsteni ri 998). Similarl y, the fastest radio m illisecond 
pulsar, PSR J1748-2446ad (|Hessels et al.l 120061). and the 
faste st nonpulsating LMXB, 4U 1608—52 (Har tman et al.l 
2003). spin at frequencies of 716 Hz and 619 Hz respec- 
tively. The gap below the breakup frequency is explained if 
the star is deformed by one part in ~ 10 8 , such that grav- 
itational radiation balances the accretion torque at hecto- 
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hertz frequencies (|Bildstenl |l998l ) . Several physical mech- 
anisms can produce the requisite deformation: ma gnetic 
mountains (IPavne fe Melatosfl2004l; iMelatos fc Pavnd l 2005l : 
iPavne fc Melatosl 120061 ; IVigelius fc Melatosl 120081 ), thermo- 
composi tional mountains c aused by electron capture gra- 
dients (lUshomirsky et all l2000l ) , toroid al internal mag- 
netic fields (ICutlerl I2OO2T), and r-mode s (|Anderssonl 1 19981 ; 
lOwen et al.ll 19981 ; iNavvar fc Owenll2006T l . AMSPs are there- 
fore promising targets for ground-based, long-baseline inter- 
ferometers like the Laser Interferometer Gravitational- Wave 
Observatory (LIGO). An AMSP at a distance of 1 kpc, spin- 
ning at 0.4 kHz with ellipticity e = 10 -8 , generates a wave 
strain h ~ 10 -27 . By comparison, initial LIGO's sensitivity 
thre shold in the 0-1-0-4 kHz band is ~ 10" 26 during the S4 
run l| Abbott et al.ll2007l ). Advanced LIGO will get down to 
h ~ 10~ 27 in the same band, and narrowband tunability will 
increase its sensitivity to AMSPs further, as f2» is known a 
priori from X-ray timing. 

An AMSP with ellipticity e ~ 10" 8 is expected to 
precess with a period of hours to days. Magnetic moun- 
tains, for example, are built around the magnetic axis, which 
is misaligned i n general with the rota tion axis in objects 
which pulsate (|Pavne fc Melatosl [2006). More generally, a 
mass quadrupole of any provenance should be kicked out 
of alignment continuously by stochastic accretion torques 
l|jones fc Anderssonl |2002| ). Hence AMSPs are promising 
observational candidates for observing short-period preces- 
sion. Until now, however, precession has been difficult to 
detect in neutron stars. Only one source, the radio pul- 
sar PSR B1 828— 11, precesses unambiguously, with period 
P p — 250 d I St airs e t al. 2000). Oscillatory trends in pulse 
arrival times, with periods of several days, have also been 
reported tentatively in a few other objec ts (|Melatos! |2000| ; 
iHobbs et al1l2006l ; IPavne fc Melatosll2006h . but the physical 
cause is unclear. 

Free precession consists of a fast wobble about the an- 
gular momentum vector J, at approximately the pulsar spin 
period P* = 2ir/Q* and a slow retrograde rotation about the 
symmetry axis, with period P p = 2n/Q v , which modulates 
the pulse shape and arrival times ijZimmermann fc Szed enits 
19791: lAlpar fc Pines! Il985l ; IJones fc Anderssonl l200ll . |2002| ; 
Linkl r2003). The precession frequency f2 p depends on the el- 
lipticity, e, and the tilt angle 9 (between the symmetry axis 
and J), with 

e cos 9 « ilp/il*. 



(1) 



The amplitude ratio of the gravitational wave sig- 
nal at the spin frequency and its second harmoni c 
ijZimmermann fc Szedenitsl 1 19791 ; Ijaranowski et all Il998h . 
and in the + and x polarizations, provides independent in- 
formation on e, 8, the orientation of J, and the emission 
pattern on the surface of the star. Narrowband tunability 
facilitates extraction of this information. 

In this paper, we compute theoretically the X-ray sig- 
nal from a precessing pulsar for a range of orientations and 
compare three quantities from each pulse profile to the data: 
the mean flux of the profile, the zero-to-peak pulse ampli- 
tude, and the pulse-phase residuals. We search for the signa- 
ture of precession in X-ray timing data from one particular 
AMSP , XTE J1814-338. An analogous search was carried 
out bv lAkgun et~al] i|2006l ) for the radio pulsar PSR B1828- 
11, who modelled the period residuals and pulse shapes an- 



alytically taking into account precession effects (biaxial and 
triaxial) as well as the contribution from the magnetic spin- 
down torque. The authors performed searches over a range of 
beam locations, degrees of triaxiality, tilt angles and angle- 
dependent spin-down torques, finding a wide range of pa- 
rameters which match the data. Thus, they were unable to 
constrain the shape of the star but did find that the angle- 
dependent spin-down torque contributes to the period resid- 
uals. Their method differs from ours in that, instead of fit- 
ting the shape of the residuals and comparing for each set of 
parameters, they determined the validity of a configuration 
by calculating Bayesian probability distribution functions 
for the parameters under certain constraints. 

The paper is structured as follows. Section [2] describes 
the precession model and its implementation. Sections and 
[4] characterize the predicted X-ray signal for a biaxial, pre- 
cessing pulsar with one and two hotspots respectively, specif- 
ically the relative precession phases between the flux, pulse 
amplitude, and pulse-phase. Section [5] repeats the predic- 
tions for a triaxial, precessing pulsar. Section [6] describes 
the data reduction and timing analysis of XTE J1814— 338. 
We compare the measurements with the theory in Section 
[7] and derive upper limits on e, 9, and the associated gravi- 
tational wave strain in Section [S] The limit on 9 constrains 
the relative strength of the driving and damping forces in 
the system. 



2 PRECESSION MODEL 
2.1 Equations of motion 

Three Euler angles (9, (p, ip) describe the rotation of a rigid 
body with body axes (ei,e2, 63) relative to the Cartesian 
triad (i,j,k) o f an inertial observer. We define (9,<p,ip) ac- 
cording to the lLandau fc Lifshit3 l|l969T ) convention. Con- 
sider, first, the special case of a biaxial, freely precessing 
neutron star. Let e3 be the symmetry axis and take k to lie 
along the total angular momentum J, as depicted in Figure 
[1] The angle between e3 and J is 9. The total angular veloc- 
ity of the precessing system, f2, comprises two components: 
tt rotates about J at a constant angle 9, with frequency <p, 
and ei and &2 rotate about e3, with frequency fl p — ip. The 
rate of precession is controlled by the ellipticity of the star, 
e, and 9, i.e. Q p = ef2cos(# + 9). Note that J, f2, and e3 are 
coplanar, as indicated by the shading in Figure [T] For small 
angles, one has 9 ~ (Ald/Ii) sin 9 cos 9, where Ji, Ja and 73 
denote the star's principal moments of inertia, and Aid is 
defined through h = h = h - AJ d /3 and h = I + 2A/ d /3. 
Aid is positive for an oblate star, negative for a prolate star, 
and is related to the ellipticity via e = Aid /lo- 
in terms of the above notation, the six equations of 
motion are 

9 = Qi cos ip — Q.2 sin (2) 
rj) — (fii sin tf> + Q.2 cos ip)/ sin 9, (3) 
ip = SI3 — cos 9(Qi sin ip + O.2 cos ip)/ sin 9, (4) 
Jrfii = (Ja-Iajnafis, (5) 

m 2 = (Js-jooiOs, (6) 

/ 3 ^3 = (Jl-/ 2 )fil«2. (7) 

The time origin is arbitrary, so initially we can set ip(0) = 



Does the accreting 




Figure 1. Geometry of precession. Drawn are the inertial frame 
axes i,j,k; the symmetry axis, e3, which makes an angle 8 with 
the total angular momentum vector, J; the angular velocity vec- 
tor, Si, which makes an angle 9 with J; and the line-of-sight vector 
n, which makes an angle \ with J. The dashed line indicates that 
n lies in the plane containing i and k. The dashed-dotted lines 
indicate that e3 and Q are coplanar. The shaded plane containing 
e3 , Si , and J rotates about J with angular frequency <f>. 



without loss of generality. The angular velocity can be 
decomposed into CI = <pk+lj)e3. This gives tp = — e^)/o/^3- In 
biaxial precession, 8 is constant. Solving ©-(El with 8 = 
and ip(0) = yields the following expressions, which we use 
to initialise the angular velocity: 

fii(0) = 0, (8) 
O 2 (0) = fi c sin0, (9) 
fi 3 (0) = n c (cos6»-eJo// 3 ). (10) 

Without loss of generality, we fix the line-of-sight vector, 
n, to lie in the i-k plane, making an angle \ with J and 
intersecting the stellar surface intensity map at latitude 67b 
and phase <f>B- The angles (#b,0b) are defined with respect 
to the moving body frame, as opposed to the Euler angles, 
which are defined with respect to the inertial frame. We then 
compute the observed intensity 7 as a function of time from a 
specific su rface intensity map F(8b, 4>b), i.e. I = F(8b,4>b)- 
Note that I Jones fc Anderssonl l|200ll ) assumed ^ x ^ 7r /2 
but did not investigate the dependence of I on \. In this 
paper, we show that the x dependence is significant. 

The above initialization gives four searchable parame- 
ters: 8, the initial azimuth 0(0), the inclination angle \i an d 
the latitude a of the hotspot(s), defined in Section \2. 21 

2.2 Synthetic pulse profiles 

Studies of the harmonic ratio and modulation amplutude 
of type I burst oscillations from six LM XBs point to the 
existence of a single, hem ispheric hotspot l|Muno et al ] |2002l : 
iPavne fc Melatoa 120061 ). We thus perform simulations for 
one and two hotspots, corresponding to 
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B (degrees) 

Figure 2. Top: Surface intensity map of the single hotspot con- 
figuration (equation llll l. Bottom: Surface intensity map of the 
double hotspot configuration (equation 1121 . 



F(8 B , 0b) = sin(<9 B + a) sin(0 B ) + DC (11) 
and 

F(8 b ,(/>b) = sin(6» B + a) sin 2 (0 B ) + DC (12) 

respectively, where DC represents a constant offset, and a 
is the latitude of the hotspot 's centre, defined to be zero at 
the equator and ±90° at the poles. For a — 0, this corre- 
sponds to a surface intensity map containing a bright spot on 
the ei-e3 plane, centred on the equator at body coordinates 
(#b,0_b) =(90°, 90°), and a dark or bright spot diametri- 
cally opposite at body coordinates (90°, 270°) as per Figure 
[5] For a 7^ 0, the body coordinates in (|lip and (|12|l are 
rotated by a in the ei-e3 plane. 

The light curves are generated by evaluating the in- 
tensity of the point on the surface map which is inter- 
cepted by the line-of-sight at each time step. These syn- 
thesised light curves are analysed in exactly the same way 
as the observational data. The profile of each pulse in the 
time series is fitted with a first and second harmonic, i.e. 
a + 6sin(27T7 + c) +dsin(47T7 + e), where a is the mean flux, 7 
denotes the pulse phase (0^7^ 1), b and d are amplitudes, 
and c and e are phase offsets. The fractional root-mean- 
square (RMS) amplitudes of the first and second harmonic 
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Figure 3. Sample simulation output for a biaxial star with a 
single hotspot. Top: mean flux, a. Center: fractional RMS, b/\/2a. 
Bottom: pulse-phase residuals, 70. The precession period is 1500 
time units. 



are b/yla and d/y/2a respectively. We define two reference 
pulse phases in each profile as the phases which maximise 
the first and second harmonic components. The pulse-phase 
residuals are the differences between the predicted and ob- 
served maxima, i.e. 70 = 1/4 — c/2n and 71 = 1/8 — e/4-7r. In 
Section [6751 we show that the second harmonic does not con- 
tribute to the overall candidate precession signal at a level 
that can be detected in the X-ray timing data presently 
available for XTE J1814-338. The quantities d/V2a and 71 
are therefore neglected in the following analysis. 

We fold the time series a, b/^/2a and 70 at the theo- 
retical precession period. These three quantities are mod- 
ulated at the precession frequency due to the motion of 
the hotspot (s) relative to the observer (see Figure [3] for an 
example of the output). The average trend in each quan- 
tity is fitted with a sinusoid, viz. A n + B n sin(27rF + C„), 
where F is the precession phase, not to be confused with 
the pulse phase, and n refers to mean flux, fractional RMS 
or pulse-phase residuals. This yields three main quantities 
of observational interest: the relative precession phase of 
a and b/y>2a (denoted by AF flux _ rms = |C flux - C rms |), 
the relative precession phase 70 and b/V2a (denoted by 
Ar phasc _ rms = ICphasc-Crmsl), and the zero-to-peak ampli- 
tude of the pulse-phase residuals (-Bphasc)- Our simulations 
run for (on average) three precession periods. 



or antiphase, i.e. Arfl ux _ rms = or n. -Bphasc increases with 
the tilt angle 9. 

The effect of varying the four searchable parameters is 
now discussed in detail. 



3.1 Tilt angle, 6 

We cons ider small til t angle s (9 00°), a s in PSR 1828-11 
[9 < 3°; IStairs et al] (|200Ch : iLinkl (|2003h ] and other astro- 
physical bodies. We verify that the precession period in- 
creases with 9 according to {TJ. We also find that Bphaae in- 
creases linearly with 9 in the small tilt angle regime 9 ^ 10° , 
ranging from 0.003 at 9 = 1° to 0.028 at 9 = 10°. There is 
also a dependence of -Bphas e on a (see bottom panel of Fig- 
ure O . [jones^&^ndersso^ (|200lf ) predicted analytically, for 
a radio pulsar with a directed beam rather than a hotspot, 
that the pulse-phase residuals vary sinusoidally on the pre- 
cession time-scale with amplitude B p h a sc oc 9/ tan(7r/2 — a'), 
where a' in their case is the latitude of the beam. This for- 
mula agrees with our results. 

Varying 9 does not affect Ar flux _ rms or AF phasc _ rms . 
Arfl ux _ rms undergoes a phase shift of n around \ ~ 90° 
which is explained below. Contour maps of AFfl ux _ rms , 
Ar p haso-rms and -B p hasc as functions of x an< i $ are plot- 
ted in Figure [4] 



3.2 Hotspot latitude, a, and inclination, \ 

Figure \S\ displays contour maps of Ar p hasc-rms, AF flux _ rms , 
and Bphasc as function of % an d a f° r 9 = 1°. In the range 
— 15° < a < 15°, Trms undergoes a n phase reversal be- 
tween certain values of x> which shows up in Ar p hasc-rms 
and AFfl ux _ rms . As \ is increased, a strong second harmonic 
gradually appears in Trms over an interval of w 10° . As x in- 
creases further, the harmonic disappears and Trms is shifts 
by 7r. The value of \ a t which this happens depends on a 
(see top and middle panels of Figure [5}. For example, at 
a = —10°, the harmonic appears for 56° < x < 64°. 

The value of AFfl ux _ rms also depends on which hemi- 
sphere (north or south) the hotspot is in, relative to n. If 
both are in the same hemisphere, then AFfl ux _ rms ~ n. If 
they are in different hemispheres, then AFfl ux _ rms ~ 0. This 
change occurs because the mean flux profile 'flips'. For exam- 
ple, if x — 30°, the fractional RMS is the same whether the 
hotspot is in the north or south. However, if the hotspot is 
in the north (south), n starts off closer to (further from) the 
hotspot, and the flux dims (brightens) as the star precesses. 
The phase shifts of Frms and the mean flux are explained 
geometrically in more detail in Section \3. 41 

-Bphasc is larger when the hotspot is close to the poles 
than when the hotspot is near the equator, but hardly varies 
with x- When x = 0, the observer does not see any evidence 
of precession at all in the pulse-phase residuals, as e3 re- 
mains equidistant from n at all times. 



3 BIAXIAL STAR: SINGLE HOTSPOT 

In most configurations involving a single hotspot, the pre- 
cession phases of b/^/2a and 70 differ by Ar p hase-rms ~ t/2. 
The precession phases of 6/v2a and a are either in phase 



3.3 Initial longitude, 0(0) 

0(0) determines the initial latitude where n intersects the 
surface. When n is within 1° of J, i.e. x < 1° or X > 179°, 
as one various 0(0) from to 2n, Ar p hasc-rms oscillates 
sinusoidally around 7r/2, AFfl ux _ rms oscillates sinusoidally 
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around or n (depending on the hemisphere of the hotspot), 
and Bphaae peaks at <j>(0) — it or 2ir. At these extremes, n 
traverses a very small area on the star's surface during each 
spin period, i.e. B rms , -Bphasc <C 1. Hence any variation due 
to <^>(0), insignificant for other values of Xi now dominates. 

At x — 90°, AFfl ux _ rms varies gradually from to n. 
This is due to the flip in the mean flux profile mentioned in 
Section 13.21 



3.4 Geometry of the phase shifts 

In this section, we explain geometrically why the relative 
precession phase of maximum rot at ion- averaged intensity 
abruptly changes whenever the line of sight crosses the equa- 
tor, causing AFfl ux _ rms to jump by 7r rad. To understand this 
counterintuitive effect, consider Figure [6] which shows the 
coplanar vectors J, e3 (which makes an angle 8 with J), and 
tt (which makes an angle 9 with J). In reality, 9 is small, but 
we enlarge it artificially for illustrative purposes. For a = 0, 
the hotspot is at the equator (relative to 63), and the sur- 
face intensity changes from brighter to darker than average 
at the dot-dashed line. Figure|6][a) is a snapshot taken of the 
pulsar at the start of a precession cycle, while Figure [6jb) is 
taken half a precession cycle later, after a time n/ip elapses. 
Both snapshots are taken at the same arbitrary pulse-phase, 
when e3 and fl line up as shown. 

Now consider two observers, Xi just north of the equator 
and X2 just south of the equator. The dashed lines show the 
loci of points where the lines-of-sight intersect the surface 
during one spin (not precession) period. The lines are tilted 
with respect to the horizontal at angle 9, perpendicular to $7, 
because the pulsar rotates instantaneously about ft during 
one spin period. At the start of the precession cycle [panel 
(a)], observer xi traces a path which passes through more of 
the dark hemisphere than the bright, while X2 traces a path 
through more of the bright hemisphere. Half a precession 
cycle later [panel (b)], the opposite happens. This means 
that the mean flux seen by xi increases from minimum to 
maximum from (a) to (b), whereas the mean flux seen by 
X2 decreases from maximum to minimum. This 'flip' causes 
the 7T phase shift in Arfl ux _ rms . 

The phase reversal in Trms at certain inclination angles 
when the hotspot is in the range —15° < a < 15° can also 
be understood with the help of Figure HJ] For a — 0, the 
line dividing the bright and dark hemispheres is at the same 
position in (a) and (b). Hence, the pulse amplitude b and 
times-of-arrivals (TOAs) seen by either observer Xl or X2 
is the same at epoch (a) and epoch (b), even though the 
mean flux a changes from epoch (a) to epoch (b). As we are 
measuring 6/v2a, the phase reversal in the mean flux causes 
a phase reversal in Frms at x — 90°. This effect occurs at 
larger (smaller) x f° r a > 0°(a < 0°), and does not occur 
when a > 15°. This is because, as the hotspot moves away 
from a = 0°, the pulse amplitudes b seen by xi and X2 
become increasingly different at epochs (a) and (b), so the 
effect of a on the fractional RMS is reduced. 

We also observe a strong harmonic component in the 
phase residuals at a — 0° as the pulse TOAs are identical 
in epochs (a) and (b). 
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Figure 4. Relative precession phase of the RMS and pulse-phase 



residuals, AR 



phasc- 



(top), RMS and mean flux, Ar fl 



ux — rms 



(middle), and amplitude of the folded pulse-phase residuals, 
Sphasc (bottom), versus tilt angle 8 and inclination angle x> both 
measured in degrees for a biaxial pulsar with one hotspot. Param- 
eters: 0(„) = 0°, a = 45° for 1° < 6 < 10°, and 5° < \ < 175°. 
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Figure 6. Schematic illustrating the rotation of one hotspot, from (a) the start of one precession cycle, to (b) midway through the 
precession cycle. The bright and dark hemispheres are marked. \\ is the position of an observer just north of the equator (x < 90°). \2 
is the position of an observer just south of the equator (x > 90°). The dashed lines indicate the points on the surface that the observers 
see each spin period. These lines are perpendicular to fl and hence tilted by the angle (artificially enhanced for clarity) with respect 
to the horizontal. The hotspot centres are indicated by the black dots. 



4 BIAXIAL STAR: TWO HOTSPOTS 

To simulate two diametrically opposed hotspots, the inten- 
sity map is changed to equation (|12[) . The pulse profiles now 
contain an increased harmonic component in certain config- 
urations. The parameter study in Section [3] is repeated. 

4.1 Tilt angle, 9 

Figure [7J shows contour maps of the three quantities as a 
function of 9 and \ f° r hotspots at 45° and 225° . As with one 
hotspot, only B p hase is affected by changing 9; it ranges from 
0.002 to 0.075 (see bottom panel of Figure [JJ. dB phllsc /d9 is 
smaller in the range 45° < x ^ 135°. 

4.2 Hotspot latitude, a and inclination, \ 

The positions of the hotspots relative to n affect 
Ar p hasG-rms and Arflux-rms- Figure [8] shows contour maps 
of the two quantities and Bphasc as functions of \ an< i 
a for 9 — 1°. We find Ar p hase-rms ~ 7t/2 in the re- 
gions where \ is north of a hotspot in the northern hemi- 
sphere, or south of one in the southern hemisphere, i.e. if 
a = 45°, 225° and x < 45° or \ > 135°. Elsewhere, we 
find Ar p hasc-rms ~ 37r/2. As for one hotspot, the fractional 
RMS reverses phase when n lies within ~ 5° of the hotspots, 
causing a second harmonic component to develop in the frac- 
tional RMS. 

As x approaches the equator, the harmonic component 
in the pulse profiles increases and peaks at x — 90°. The 
horizontal band across the contour map at x = 90° is caused 
by the pulse-phase residuals shifting by n (see top panel of 
Figure [5}. At a = 0°, the two hotspots are within 9 of n. For 
small 9, the observer sees the fractional RMS, mean flux and 
pulse-phase residuals rise and fall twice as the star precesses 
about e3 (at rate ip)- This halves the apparent precession 
period. 

Arfl ux _ rms varies from to n in a similar fashion. The 
fractional RMS and mean flux are generally out of phase 
when n is north of a hotspot in the northern hemisphere, 
south of a hotspot in the southern hemisphere, and vice 



versa. However, the fractional RMS and mean flux undergo 
phase shifts at different points (x> a )- The fractional RMS 
'flips' at the points indicated by the diamond pattern in 
the top panel of Figure [8] whereas the mean flux 'flips' for 
65° < x IS 115°, indicated by the horizontal band in the 
middle panel of Figure [8] The origin of these phase shifts is 
explained geometrically in Section H. 41 

-B p haso is largest when both x an d « are close to the 
poles. 

4.3 Initial longitude, 0(0) 

Longitudinal dependences are weak. For example, 
Ar p haso-rms varies by < 0.02% as 0(0) goes from 0° 
to 330° for x = 70°, as is shown in Figure M Near x = 90°, 
Ar ph aso-rms and AF flux _ rms vary more strongly with 
0(0) due to phase shifts in the pulse-phase residuals and 
fractional RMS at various configurations. However, these 
shifts do not translate into true changes in the observed 
pulse profile, because the second harmonic component, 
which we do not consider, dominates in this region and 
behaves differently. The mean flux profile is not affected. 
Arfl ux _ rms cycles between ~ 0.9 and 27r as 0(0) varies from 
to 2-7T rad, while Ar p h aS c-rms cycles between 7r/2 and 
3tt/2. 

4.4 Geometry of the phase shifts 

Figure [TU] shows a pulsar with two hotspots at (a) the be- 
ginning and (b) midway through a precession cycle. The 
hotspots are located at approximately 45° , 225° to match 
Figures [7J and [5] The shaded band indicates the region of 
the surface which is darker than average; the darkest points 
lie along the dot-dashed line. The band changes orientation 
from (a) to (b) as the star precesses rigidly about e3 (after 
a time ir/ip elapses). 

To explain the flip in the fractional RMS at x — 45° (or 
equivalently 135°), we compare the observers at xi (north of 
45°) and X2 (south of 45°). In (a) and (b), Xi traces similar 
paths close to the darkest and brightest regions respectively, 
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Figure 5. Relative precession phase of the RMS and pulse-phase 
residuals, Ar p h a sc-rms (top), RMS and mean flux, Arfl ux _ rms 
(middle), and amplitude of the folded pulse-phase residuals, 
^phasc (bottom), versus inclination angle \ an d hotspot latitude 
o, both measured in degrees for a biaxial pulsar with one hotspot. 
Parameters: 9 = 1°, (o) = 0, 5° < a < 175°, and 5° < \ =S 175°. 



Figure 7. Relative precession phase of the RMS and pulse-phase 
residuals, Ar p i, asc _ rlns (top), RMS and mean flux, Ar Hux _ rms 
(middle), and amplitude of the folded pulse-phase residuals, 
Sphasc (bottom), versus inclination angle \ and tilt angle 8, both 
measured in degrees, for a biaxial pulsar with two hotspots. Pa- 
rameters: 0(0) = 0, a = 45°, 1° Sj 9 ^ 10° and 5° ^ \ < 175°. 
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Figure 8. Relative precession phase of the RMS and pulse-phase 



residuals, Ar 



phase — rms 



(top), RMS and mean flux, Ar fl 



ux — rms 



(middle), and amplitude of the folded pulse-phase residuals, 
^phasc (bottom), versus inclination angle \ an< J hotspot lati- 
tude o, both measured in degrees, for a biaxial pulsar with two 
hotspots. Parameters: 6 = 1°, 0(0) = 0, 5° < a ^ 175° and 
5° sS X 175°. 



Figure 9. Relative precession phase of the RMS and pulse-phase 
residuals, Ar p h asc _ rms (top), RMS and mean flux, Arfl ux _ rms 
(middle), and amplitude of the folded pulse-phase residuals, 
Sphasc (bottom), versus inclination angle \ an d initial phase 
</>(0), both measured in degrees, for a biaxial pulsar with two 
hotspots. Parameters: 9 = 1°, a = 45°, 0° < (f> < 330° and 
0.1° < x < 179.9°. 
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leading to similar pulse amplitudes. However, the mean flux 
is lower in (a) than in (b), so the fractional RMS is a max- 
imum at (a) and a minimum at (b). Similarly, X2 traces 
a path with a greater pulse amplitude in (a) than in (b), 
but because these are normalised by the mean flux, which is 
greater at (a) than at (b), the fractional RMS is a minimum 
at (a) and a maximum at (b). 

In order to explain the flip in mean flux between 65° < 
X < 115°, we compare X2 (south of 65°) and X3 (south of 
115°). In Figure [TUT a). both X2 and X3 trace paths of similar 
brightness over one spin period. However in Figure fTUT b). X2 
traces a path in the darker band, whereas X3 traces a path 
closer to the bright centre of the hotspot. This confirms that 
the mean flux profiles seen by X2 and X3 are T out of phase, 
as our simulations show. 



6 X-RAY TIMING ANALYSIS 

In Sections we show that precession modulates the 

mean intensity and arrival times of pulses from one or two 
hotspots. In this section, we search for such modulation in 
X-ray timing observations of AMSPs made over recent years 
by the Rossi X-ray Timing Explorer (RXTE). We consider 
three sources observed by RXTE: SAX J1808.4-3658, XTE 
J1814-338, and HETE J1900. 1-2455. XTE 1900.1-2455 
is unsuitable as it displays peculiar behaviour, including 
persistent X-ray emission and i ntermittent pulsation s even 
during periods of low accretion (|Gallowav et al.ll2007l ). SAX 
J1808.4— 3658 is also not ideal as the fractional RMS of the 
first harmonic changed erratically over the last four out- 
bursts, and a previous analysis of the pulse-phase residu- 
als di d not reveal any per iodicities consistent with preces- 
sion l|Hartman et al.l 20071 ) . The most promising candidate 
is XTE J1814— 338, whose data cover a 66-day outburst in 
which modulations in the mean flux, RMS, and pulse-phase 
residuals are visible by eye (see below). 



5 TRIAXIAL STAR 



A triaxial star is described by two separate ellipticity pa- 
rameters, ei and 62, d efined by ei = [2(73 — JiV/il 1 ^ 2 an d 
£2 = [2{h~h)/h] 1/2 (|Zimmermann fc Szedenitsll 19791 ). We 
pick Ji to have a reference value of Jo = 0.4M*P 2 and in- 
vestigate the effect of varying 62/61 from 0.3 to 0.7. Within 
this range, we find that there are no changes in the be- 
haviour of Ar phase _ rms , Ar flux _ rms or Bphasc- The only ef- 
fect is to change the spin period, i.e. P0.3 ~ 1.0022Po.5, 
Po.7 ~ 0.9994Po.5 where P0.3 is the spin period for 62/61 = 

0.3 and so on. 

The motion is discussed in detail bv lLandau fc Life hitz 

(1969). When J 2 is only slightly larger than 2Eh, where E 
is the total energy, the ei axis rotates around J in an ellipse 
whose size increases as J increases. As J 2 approaches 2ET>,, 
e3 rotates around J in an ellipse. Hence 9 is not constant: 
its mean, (8), increases with J. 

We consider small 9 in astrophysical problems, where 
9 now denotes the tilt angle at the start of the simulation. 
The angular velocity vector fl rotates periodically around 
the body axes, with period T ~ 2ir/(eQ) (the exact value 
is given by a complete elliptic integral of the first kind). 
However, ft does not return to its original position with 
respect to the inertial axes after one cycle. For a biaxial 
body, the spin frequency is simply <j>- For a triaxial body, we 
have <f>(t) = 4>i{t)+4>2{t), where <j>i(t) has period T and <fe(£) 
has period T' [incommensurable with T; |Landau fc Lifshita 
1 19691 ; IZimmermann fc Szedenitelll979l )] . 

We repeat the parameter searches in Sections [3] and [4] 
and find similar results for one and two hotspots. In some 
configurations, the precession phase profiles of Ar p h asc _ rms , 
Arfl ux „ rms and -Bphase become non sinusoidal or contain 
strong harmonics. For consistency, we fit these profiles with 
the same format as before (see Section 12. 2\ , recovering a 
modulation similar to Figures [TH9] 

The only noticeable difference is in the case of two 
hotspots. At x = 90°, between a « ±30° to ±50°, 
AFfl ux _ rms increases from 3.14 to 3.19. There is also in- 
creased variation with (f)(0) (< 5% for Ar phasc _ rms and 
< 7.5% for AF fl 

ux-rras ) • 



6.1 XTE J1814 338 

XTE J1814-338 is the fifth A MSP to be discovered wit h 
a spin frequency of 314.4 Hz (|Markwardt fc Swankl I2003T ) . 
Between MJD 52796 and 52834, the object experienced 
an ouburst during which 27 thermonuclear (type I) X- 
ray bursts were observed. This is the longest interval over 
which pulsations have been detected consistently. The data 
from this outb urst were analysed previ ously, but with differ- 
ent em phases. IWatts et all l|2005l ) and I Watts fc Strohmaverl 

(2006) reported on the variability and energy dependence 
of these bursts, finding that the burst fractional amplitude 
(defined in the above papers) is constant d uring a burst 
and d ecreases with increasing photon energy. iPapitto et al.l 

(2007) presented a timing analysis and refined orbital pa- 
rameters previously published. They noted the modulation 
in the pulse-phase residuals and attributed it to movement 
of the accretion hotspot as the accretion rate varies. 



6.2 Observations 

The source was observed b y the RXTE Prop ortional 
Counter Array [RXTE PCA; [jahoda et alj (|l996T )]. which 
consists of five proportional counter units (PCUs). During 
the course of an observation, different numbers of PCUs are 
turned on at different epochs, even within one data set. 
In order to accurately determine the background rate for 
each data set, the contribution from each PCU at all times 
must be tracked and accounted for. Data were collected in 
Event Mode over 64 energy channels (2-60 keV) with 125 /j,s 
time resolution. The data comprise 91 internally contiguous 
blocks lasting from 2 to 30 ks and span a total of 66 days. 
The X-ray flux was measured by fitting a phenomenological 
model consisting of blackbody and power-law components, 
each attenuated by neutral absorption, to spectra extracted 
from Standard-2 mode data in the range 2.5-25 keV. The 
Standard-2 mode has 129 energy channels and 16 s timing 
resolution. 
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Figure 10. Schematic illustrating the rotation of two hotspots, from (a) the start of one precession cycle, to (b) midway through the 
precession cycle. The shaded band indicates the region on the surface that is darger than average, for a R- 45°. xi i s the position of 
an observer just north of a hotspot (\ < 45°), \2 is the position of an observer south of a hotspot (\ > 65°). X3 i s the position of an 
observer south of 115°, where the mean flux profile is seen to flip. The dashed lines indicates the points on the surface that the observers 
see each spin period. These lines are perpendicular to f2 and hence tilted by the angle 9 (artificially enhanced for clarity) with respect 
to the horizontal. The hotspot centres are indicated by the black dots. 



6.3 Timing analysis 

The data is processed using LHEASOFTl] version 5.3 (2003 
November 17). We correct the photon times of arrival 
(TOAs) to the solar system barycentre. The X-ray flux is 
background subtracted using the RXTE /PC A Mission- Long 
Bright Source background mode0. A separate response ma- 
trix is calculated for each observation to account for drift in 
the PCA gains. We remove all type I bursts from the data, 
subtracting all photons from 15 seconds before to 200 sec- 
onds after the burst peak; all the bursts have r ise times of 
1-8 s and last from 100-200 s l| Watts et al.ll2005h . 

The X-ray flux is plotted versus time in Figure 1111 It 
rises over the first five days from 3.5 x 10 -10 to 4.4 x 10~ 10 
erg cm" s~ . It remains between 4.4 x 10~ 10 and 5.1 x 10" 10 
erg cm -2 s -1 for the next 30 days, exhibiting strong modula- 
tions. It then drops sharply over the next two days and falls 
below the sensitivity threshold (0.2 xl0~ 10 ergcm~ 2 s _1 ) 
on MJD 52834. We analyse the first 37 days of data only 
(between the dashed vertical lines in Figure [TTjl . 

The published orbital parameters at th e time of writing 
l|Markwardt et al"1l2003l ; IChung et al.ll2006l) are incomplete: 
neither reference quoted the epoch of mean phase (tgo). To 
extract this quantity, we select a long barycentre-corrected 
data span (observation ID 80418-01-03-00, lasting 30.024 ks) 
and correct for the satellite orbit using a trial tgo value. Any 
error in the trial tgo value, or any other orbital parame- 
ter, can be extracted by comparing the actual and expected 
TOAs . The TOA residuals obey equation (3) in lDeeter et al.l 
|l98ll ). 

The fifth term on the right-hand side of the latter equa- 
tion shows that the residuals in £90 produce a sinusoidal 
variation in the TOAs, which we attempt to minimize. We 
also calculate corrections to the spin period by subtracting a 
linear trend from the TOAs, and the spin period derivatives 



1 http:/ /heasarc. gsfc.nasa.gov/lheasoft 

2 http:/ /heasarc. gsfc.nasa.gov/docs/xte/pca_news. html 



Table 1. Orbital parameters for XTE J1814-338 

Barycentric spin period (s) 0.003181105669954(4) 

Spin frequency derivative (Hzs - - 1 ) —7.2(3) X 10 -14 

Projected semimajor axis (It s) 0.390626(2) a 

Epoch of tt/2 mean phase (MJD) 52808.8975258(4) 

Orbital period (s) 15388.7238(2) a 



a IChung et all <|2006T ) 



by sub tracting a quadratic trend [third term in lDeeter et al.l 
(|l981l )]. We do not attempt to correct for other effects as 
these three terms dominate. The revised orbital parameters 
are quoted in Table [1] 



6.4 Pulse folding 

We fold the light curves for each contiguous observing seg- 
ment at the measured spin period for that segment, taking 
into account the first-order spin frequency derivative. 

The folded pulse profiles are then fitted with a sinusoid 
comprising first and second harmonic components, as with 
the simulations above, viz. a-r-bsin(27T7 + c) +dsin(47r7 + e), 
where ^ 7 ^ 1 denotes the pulse-phase. Fitting is done 
using the Levenberg-Marquardt nonlinear least-squares al- 
gorithm. Again, from the fitted parameters, we measure the 
fractional RMS of the first and second harmonic (fo/\/2a 
and d/y/2a respectively) as well as the pulse-phase residu- 
als of the first and second harmonic (70 = 0.25 — c/27r and 
71 = 0.125 — d/47r respectively). 

Uncertainties in the fitted paramet ers are determine d 
using the constant \ 2 boundary method (|Press et al.lll986t ); 
i.e. b, c, d and e are iterated separately until x 2 increases by 
unity relative to its minimum. We find that the constant 
X 2 uncertainties are ~ 1.5 times the raw standard deviation 
from the least-squares fit, a^ t . Henceforth, to simplify the 
extensive analysis, we use 1.5o"fit to qualify the uncertainties. 
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Figure 11. Total flux (in units of 10~ 9 erg/cm 2 /s) versus time 
(Modified Julian Date). Top: The entire 66 day outburst starting 
at MJD 52796. Bottom: The first 37 days. The flux is measured 
with the RXTE PCA in the energy band 2.5-25 keV. It is plotted 
with 1-cr error bars. Overplotted is a quadratic polynomial which 
is later subtracted in order to extract the amplitude and period 
of the ~ 12-day modulation (the 'bumps' in the graph). 

6.5 First harmonic 

In order to clarify whether to include the second harmonic 
of the observed pulse profile in any further analysis, we take 
the fitted parameters of the first harmonic to be the 'true' 
parameters and investigate the contribution of the second 
harmonic to the total signal. In the first 37 days of data, the 
first and second harmonic pulse-phase residuals 70 and 71 
are similar. There is a 1.04% difference in their gradients, 
and a linear trend between the two can be fitted with a slope 
of unity lying within the 4a limit. As for the fractional RMS, 
the first harmonic b/y/2a shows a 4.19 x 10 -4 day -1 increase 
over the first 37 days, whereas the second harmonic d/y2a 
decreases by 3.94 x 10 -4 day -1 . 

We calculate the Lomb-Scargle periodogram l|Lombl 
ll976J : |Press et al.lll986l ) for 70, 71, b/V2a and d/V^a. The 
periodogram is discussed fully in Section 16.61 For now, we 
merely note that a significant (> 4.5a) 12.2-day periodic sig- 
nal is found in 70 and b/V^a. In the second harmonic, this 
signal is present in 71 at a 3.5a level, but absent in d/v2o. 



In order to test the quality of the data, we fit 12.2-day sine 
waves to b/y/2a and d/y/2a. The 3a upper limit on the frac- 
tional amplitude of the fit in d/\/2a is 7.8%, whereas the 
best fit fractional amplitude in b/y/2a is 6.9%. We cannot 
therefore rule out the possibility that there is a hidden signal 
in the second harmonic fractional RMS. 

For this reason, and since the first harmonic component 
of the fractional RMS dominates (b/V^a « 0.103, d/y/2a « 
0.029), we exclude the second harmonic fractional RMS from 
further consideration. We can also exclude 71 as it offers 
no additional information about the candidate precession 
signal. Nonetheless, in another object with a cleaner signal, 
or with better data, the pulse-phase residuals can yield extra 
information in principl e, e.g. about the deta iled form of the 
surf ace intensity map jHartman et alj l2007h . 

lHartman et al.l (|2007f ) discussed how 'red noise' (long- 
time-scale correlations) affect each harmonic's fractional 
RMS and pulse-phase residuals, causing them to vary inde- 
pendently. The authors use a common phase residual which 
is a weighted combination of 70 and 71 , to minimise the in- 
trinsic variation. We do not employ this technique as the 
same modulation appears in 70 and 71, and it is not de- 
tectable in the second harmonic fractional RMS. 

6.6 Searching for precession 

The three remaining quantities of interest are the total flux, 
fractional RMS and pulse-phase residuals of the first har- 
monic. As any longer period variation can mask shorter pe- 
riods, a quadratic trend in the flux and a linear trend in 
the fractional RMS are subtracted, leaving the time series 
in Figure 1121 To search for periodicities in these quantities, 
we construct a Lomb-Scargle periodogram, which calculates 
the significance of periodiciti es in unevenly sampled data 
|Lomblll976l ; iPress et all 119861 ) . The periodogram is plotted 
in Figure ITUl The peak Lomb power is 18.1 for 70, 16.2 for 
6/v2a, and 16.4 for a. These Lomb powers correspond to 
a significance (the probability of a falsely detected signal) 
of order 10 -6 for 70 and 10 -5 for the other two quantities. 
The peak Lomb powers occur at periods Pfl ux = 283 ± 22 
hours, Pphase = 293 ± 21 hours, and Prms = 302 ± 23 
hours respectively (refer to Figure [13}. The uncertainties in 
these periods are obtained in Monte-Carlo fashion by adding 
quasi-random noise (a Gaussian distribution with the same 
standard deviation as the original data minus a pure sine 
wave with the respective periods Pa ux , Pphase and Prms) to 
noiseless signals with periods Pfl ux , Pphase and Prms- 

The three periods are consistent and imply a mean can- 
didate precession period of P p — 293 ± 22 hours. 

We now fold the flux, fractional RMS, and pulse-phase 
residuals at P v . The results are shown in Figure 1141 We 
describe the resulting folded time series, as with the simu- 
lations, with a sinusoid Bsin(27rr + C) (where ^ Y ^ 1 
is the precession phase) plus a DC offset. The fitted param- 
eters are listed in Table [2] From these parameters, we can 
construct three quantities which are independent of mean 
flux and determined only by the orientation of the pulsar: 
the relative precession phase offset between the flux and the 
fractional RMS, AFfl ux _ rms , the relative precession phase 
offset between the pulse-phase residuals and the fractional 
RMS, AFphaso-rms; and the amplitude of the folded pulse- 
phase residual profile, Pphase - The measured values are listed 
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Figure 12. Top: flux time series (folded on the spin period) af- 
ter subtracting the long-term quadratic trend in Figure [TT] Cen- 
ter: fractional RMS of the first harmonic component of the folded 
pulse profile after folding on the spin period. Bottom: Pulse-phase 
residuals of first harmonic. All quantities are in the energy band 
2.5—25 keV and are plotted with 1-u error bars, derived as ex- 
plained in Section [52] 



Table 2. Results of fitting the folded mean flux, fractional RMS, 
and pulse-phase residuals (folding period = 292.59 hours) with 
the expression A + Bsin(27iT + C), where T is the precession 
phase. We list the parameters relevant to the precession model, 
i.e. the amplitude (B TrnB ) and precession-phase offset (C rms ) of 
the fractional RMS, the amplitude (£? p hase) an d precession-phase 
offset (Cphase) of the pulse-phase residuals, the precession-phase 
offset of the mean flux (Cfl ux ), the relative precession phase of 
the fractional RMS and pulse-phase residuals (Ar p ha S c-rms); an d 
relative precession phase of the fractional RMS and mean flux 
(Arg ux _ rms ). Phases are expressed in radians. 
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Figure 13. Lomb periodogram for the following time series. Top: 
mean flux. Center: fractional RMS. Bottom: pulse-phase residu- 
als. Frequencies are measured in Hz. 



in Table [2] Errors are obtained using the constant \ 2 bound- 
ary method described above. 

In addition to the quantities in Tables [5] there is one ar- 
bitrary precession phase (corresponding to a choice of time 
origin) and two arbitrary amplitudes (for the mean flux and 
RMS), which depend on the unknown (and possibly chang- 
ing) DC offset flux. The latter offset in general comes from 
orientation effects and also a DC component in the surface 
intensity map. 



7 COMPARISON BETWEEN DATA AND 
SIMULATIONS 

The only configuration for which the measured values of 
Ar phasc - rm s, Ar flux _ rms and B piiaae come close to matching 
the data within experimental errors is a biaxial star with one 
hotspot. This agree s with the s i ngle-h emispheric-hotspot 
model suggested by iMuno et all l|2002l ). For 6 = 3°,x = 
179.95°, 4> = 210°, a < 0°, we find Ar phase _ rm8 = 2.6, 
Ar flux _ rms = 1.0 and B p ha SC = 0.024. We plot T flux ,r rms , 
and 70, refolded over the model's precession period, in Fig- 
ure [15] The precession phase in the plots is offset to match 
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Figure 14. Slow variation of the pulse characteristics refolded 
on the candidate precession period P v = 293 hr. Top: mean flux. 
Center: fractional RMS of the first harmonic. Bottom: pulse-phase 
residuals of the first harmonic. Measurements are represented by 
open symbols; the binned profile is represented by the histogram. 
The fitted sinusoid (with DC offset) is depicted as a solid curve. 



Figure 15. Refolded simulated time series for the best match 
configuration (9 = 3°,x = 179.95°, 4> = 210°, a = 45°). 
AFflux— rms = 1.0 and -Bphasc = 0.024 match the data, but 
^rphasc-nns = 2.6 lies slightly outside the 1 — a error bar. Top: 
mean flux. Center: fractional RMS of the first harmonic. Bottom: 
pulse-phase residuals of the first harmonic. Simulated data are 
graphed as solid curves. Best-fit sinusoid is graphed as a dashed 
curve. 



the data (cf. Figure [TJ). The exact latitude of the hotspot 
does not affect these results, as mentioned previously (e.g. if 
a < 0°, then we would get the same results for \ — 0.05°). 
The above model for XTE J1814-338 is unlikely a priori. 
Given that pulsars are oriented randomly relative to an ob- 
server, the likelihood of observing a pulsar with \ < 1° or 
X > 179° is 0.008%. In fact, 95% of the sky area covers the 
range 18° < x 162°. Secondly, for such a small inclination 
angle, n traces a small circle on the pulsar's surface during 
each rotation, resulting in a smaller (< 1%) fractional RMS 
than the ~ 12% level seen in the data. 

For configurations with 1° < x < 179°, our simulations 
fail to match the data. The relative precession phases of the 
fractional RMS, mean flux and phase residuals do not fall 
within la of the measured values. Hence, if x truly does lie 
in the above range, two scenarios are possible. 



Scenario 1 

If the 293-hr modulation in the data is a real precession sig- 
nal, our model is incomplete. For example, the intensity map 
may be more complicated in reality than equations (|ll|l or 
(|12l perhaps explaining the discrepancy between the simu- 
lations and the data in AFfl ux _ rms and AF p h asc _ rms . 

Unlike Ar flux _ rms and Ar phasc _r ma , it is possible to 
match Bphase for a relatively broad range of x- For one 
hotspot, the simulations match (0.021 ^ Bphasc ^ 0.027) 
for 6° ^ 8 ^ 10°. For a triaxial star, the match occurs at 
certain combinations of 6 and x ( see Figure [7} , with the 
most probable combination drawn from 5° ^ 6 ^ 8° and 
60° ^ x ^ 120°. Although this agreement is insufficient as 
a proof of precession without an intensity map that also re- 
produces AFfl ux _ rms and Ar phasc _ rms , it does provide some 
insight as to what the tilt angle would be in such a scenario. 

Although we did not include the type I bursts in our 
analysis, the phase residuals of the burst oscillations are 
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phase-locked with the non-accreting pulse phase residuals 
and are modulated on the same time- scale over the span 
of data that we use (| Watts et al.l [2008). This supports the 
precession model since precession of the entire pulsar would 
move the burst location(s) along with the non-accreting re- 
gions. The pulses from both areas would therefore be mod- 
ulated in the same way. 

Assuming a precession period of 293 hours and a tilt 
angle of 6°, equation |Q} implies e ~ 1CT 9 for XTE J1814- 
338. The gravitational wave strain ho at Earth from a biaxial 
rotator is given by 



ho = 



c 4 d ' 



(13) 



where G is the gravitational constant, c is the speed of 
light, Io is the star's moment of inertia and d is the dis- 
tance to the source. For XTE J181 4— 338, we have d w 
8 ± 1.6 kpc dStrohmayer et allliool ) and Rc 2 /GM > 4.2 
l)Bhattacharvva et al] 1200a ). where R is the star's radius, 
implying 10 -28 s: ho < 10 -27 . 



Scenario 2 

If the 293-hr mod ulation is not due to precession [e.g. 
iPapitto et al.l (|2007h suggested that the hotspot drifts peri- 
odically around the star] , then either the precession is heav- 
ily damped, making 9 very small, or e itself is smaller than 
expected. 

If it is rigid, the star has a nonzero ellipticity for the rea- 
sons listed in Section[T] However, in reality the star is elastic 
and probably contains a superfluid interior. Hence preces- 
sion is d amped via internal d issipation and gravitational ra- 
diation l|Cutler fc Jo nes 2001). Intern al dissipation generall y 
dominates. Based on calculations bv lAlpar fc Saula ill988h. 
the ti me-scale for damping the tilt angle iBondi fc Gold! 
l|l955h is predicted to be between 400 and 10 4 precession 
periods. This effect or a small ellipticity, or some combi- 
nation of both, lengthens the precession period beyond the 
~ 37-day observation window for ecosO < (P*/37-d), i.e. 
ecos6» < 9.9 x 10~ 10 . This implies h ^ 10 X 10 -27 cos(<9) _1 . 



to find a model configuration which matches the measured 
Ar p hasc-rms or Arflu X _ rm s, unless we choose 0° < x ^ 1° 
(or 179° ^ x ^ 180°), an a priori unlikely orientation. How- 
ever, we are able to match Bphasc for a range of tilt angles 
5° ^ ^ 10° , if we are prepared to tolerate a discrepancy 
of 50° in Ar phasc _ rms and 55° in Arfl ux _ rms between the 
data and the model. One can therefore draw two possible 
conclusions: either the star is precessing but our surface in- 
tensity map is too simplistic, or the source is not precessing. 
If we attribute the 12.2-d periodicity to precession, this im- 
plies an ellipticity of e ^ 3 x 10~ 9 , a gravitational wave 
strain ho ^ 10 -27 , and hence a signal-to-noise ratio of 10 -3 
for initial LIGO and 10 -2 for advanced LIGO (for a coher- 
ent 120-day search). On the other hand, if the precession 
is damped by internal dissipation (6 is small), or the pre- 
cession period is much longer than the 37-day data span (e 
is small), some other mechanism must cause the observed 
modulation. In this scenario, we find ecos8 ^ 9.9 x 10~ 10 
and h ^ 10 x 10~ 27 cos(6>) _1 . 

Although we face a negative result for this particu- 
lar source, this paper establishes a framework for analyz- 
ing modulations in X-ray flux from AMSPs for a range 
of geometrical configurations and surface intensity maps. 
We anticipate that the framework will be applied to other 
AMSPs in the future. Given the values of e inferred from 
the gravitational-wave stalling hypothesis dBildstenl Il998h 
and the theoretical models, e.g. of magnetic mountains 
|Pavne fc Melatosll2006l ; IVigelius fc Melatosll2008l ). it is clear 
that long-term X-ray monitoring of AMSPs (over years) is 
essential for predicting, and then searching for, their gravi- 
tational wave signal. 
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8 CONCLUSION 

By analyzing X-ray timing data from the accreting millisec- 
ond pulsar XTE J1814-338, we find a 12.2-day periodicity 
in the mean flux, fractional RMS, and pulse-phase residu- 
als of the first harmonic of the folded pulse. We measure 
two precession phase offsets relating these three quantities 
(APphaso-rms = 3.1 ± 0.2 rad and APfl ux _ rms = 0.7 ± 0.3 
rad) as well as the amplitude of the pulse-phase residuals, 
Bphasc = 0.024 ± 0.003. 

Simulations of biaxial and triaxial precessing pulsars 
with one and two hotspots are also performed for a range 
of inclination angles (0° < \ < 180°), tilt angles (6 ^ 10°), 
and hotspot latitudes (—85° ^ a ^ 85°). -Bphasc is found to 
depend on the tilt angle, while AP p hasc-rms and Arfl ux _ rms 
depend on the relative orientations of the line-of-sight and 
the hotspot (s). We find no significant dependence on the 
initial longitude at which the line-of-sight intersects the star 
except for small (< 1°) inclination angles. 

Comparing the data with the simulations, we are unable 



REFERENCES 

Abbott B., Abbott R., Adhikari R., Agresti J., Ajith P., 
Allen B., Amin R., Anderson S. B., Anderson W. G., 
Arain M., Araya M., Armandula H., Ashley M., Aston 
S., Aufmuth P., 2007, Phys. Rev. D, 76, 042001 
Akgun T., Link B., Wasserman I., 2006, MNRAS, 365, 653 
Alpar M. A., Cheng A. F., Ruderman M. A., Shaham J., 

1982, Nature, 300, 728 
Alpar M. A., Pines D., 1985, Nature, 314, 334 
Alpar M. A., Sauls J. A., 1988, ApJ, 327, 723 
Andersson N., 1998, ApJ, 502, 708 

Bhattacharyya S., Strohmayer T. E., Miller M. C, Mark- 
wardt C. B., 2005, ApJ, 619, 483 

Bildsten L., 1998, ApJL, 501, L89+ 

Bondi H., Gold T., 1955, MNRAS, 115, 41 

Chung Y.-y., Chou Y., Tsai O, Chang H., 2006, in Bulletin 
of the American Astronomical Society Vol. 38 of Bulletin 
of the American Astronomical Society, Fine Orbital Pa- 
rameters and Soft Phase Lags of Accretion Powered Mil- 
lisecond Pulsar XTE J1814-338. pp 81— I- 



Does the accreting millisecond pulsar XTE J1814—338 precess? 15 



Cook G. B., Shapiro S. L., Teukolsky S. A., 1994, ApJ, 424, 
823 

Cutler C, 2002, Phys. Rev. D, 66, 084025 
Cutler C, Jones D. I., 2001, Phys. Rev. D, 63, 024002 
Deeter J. E., Pravdo S. H., Boynton P. E., 1981, ApJ, 247, 
1003 

Eckert D., Walter R., Kretschmar P., Mas-Hesse M., 
Palumbo G. G. C, Roques J. -P., Ubertini P., Winkler 
C, 2004, The Astronomer's Telegram, 352, 1 

Galloway D., 2007, in Kang Y. W., Lee H.-W., Leung K.- 

C, Cheng K.-S., eds, The Seventh Pacific Rim Conference 
on Stellar Astrophysics Vol. 362 of Astronomical Society 
of the Pacific Conference Series, Pushing the Limit on 
Neutron Star Spin Rates, pp 105 — h 

Galloway D. K., Morgan E. H., Krauss M. I., Kaaret P., 
Chakrabarty D., 2007, ApJL, 654, L73 

Hartman J. M., Chakrabarty D., Galloway D. K., Muno 
M. P., Savov P., Mendcz M., van Straaten S., Di Salvo T., 
2003, in Bulletin of the American Astronomical Society 
Vol. 35 of Bulletin of the American Astronomical Society, 
Discovery of 619 Hz Thermonuclear Burst Oscillations in 
the Low-Mass X-Ray Binary 4U 1608-52. pp 865-+ 

Hartman J. M., Patruno A., Chakrabarty D., Kaplan D. L., 
Markwardt C. B., Morgan E. H., Ray P. S., van der Klis 
M., Wijnands R., 2007, ArXiv e-prints, 708 

Hessels J. W. T., Ransom S. M., Stairs I. H., Freire P. C. C, 
Kaspi V. M., Camilo F., 2006, Science, 311, 1901 

Hobbs G., Lyne A., Kramer M., 2006, Chinese Journal of 
Astronomy and Astrophysics Supplement, 6, 169 

Jahoda K., Swank J. H., Giles A. B., Stark M. J., 
Strohmayer T., Zhang W., Morgan E. H., 1996, in Sieg- 
mund O. H., Gummin M. A., eds, Proc. SPIE Vol. 2808, 
p. 59-70, EUV, X-Ray, and Gamma-Ray Instrumenta- 
tion for Astronomy VII, Oswald H. Siegmund; Mark A. 
Gummin; Eds. Vol. 2808 of Presented at the Society of 
Photo-Optical Instrumentation Engineers (SPIE) Confer- 
ence, In-orbit performance and calibration of the Rossi X- 
ray Timing Explorer (RXTE) Proportional Counter Array 
(PCA). pp 59-70 

Jaranowski P., Krolak A., Schutz B. F., 1998, Phys. Rev. 

D, 58, 063001 

Jones D. I., Andersson N., 2001, MNRAS, 324, 811 
Jones D. I., Andersson N., 2002, MNRAS, 331, 203 
Krimm H. A., Markwardt C. B., Deloye C. J., Romano P., 
Chakrabarty D., Campana S., Cummings J. R., Galloway 

D. K., Gehrels N., Hartman J. M., Kaaret P., Morgan 

E. H., Tueller J., 2007, ArXiv e-prints, 709 
Kulkarni A. K., Romanova M. M., 2005, ApJ, 633, 349 
Landau L. D., Lifshitz E. M., 1969, Mechanics. Course of 

Theoretical Physics, Oxford: Pergamon Press, 1969, 2nd 
ed. 

Lasota J.-P., 2001, New Astronomy Review, 45, 449 
Link B., 2003, in Bailes M., Nice D. J., Thorsett S. E., 

eds, Radio Pulsars Vol. 302 of Astronomical Society of the 

Pacific Conference Series, Precession of Isolated Neutron 

Stars, pp 241-+ 
Lomb N. R., 1976, Astrophys. Space. Sci., 39, 447 
Markwardt C. B., Strohmayer T. E., Swank J. H., 2003, 

The Astronomer's Telegram, 164, 1 
Markwardt C. B., Swank J. H., 2003, IAU Circ, 8144, 1 
Melatos A., 2000, MNRAS, 313, 217 
Melatos A., Payne D. J. B., 2005, ApJ, 623, 1044 



Morgan E., Kaaret P., Vanderspek R., 2005, The As- 
tronomer's Telegram, 523, 1 
Muno M. P., Ozel F., Chakrabarty D., 2002, ApJ, 581, 550 
Nayyar M., Owen B. J., 2006, Phys. Rev. D, 73, 084001 
Owen B. J., Lindblom L., Cutler C, Schutz B. F., Vecchio 

A., Andersson N., 1998, Phys. Rev. D, 58, 084020 
Papitto A., di Salvo T., Burderi L., Menna M. T., 

Lavagetto G., Riggio A., 2007, MNRAS, 375, 971 
Payne D. J. B., Melatos A., 2004, MNRAS, 351, 569 
Payne D. J. B., Melatos A., 2006, ApJ, 641, 471 
Press W. H., Flannery B. P., Teukolsky S. A., Vetterling 
W. T., 1986, Numerical Recipes: The Art of Scientific 
Computing, 1st edn. Cambridge University Press, Cam- 
bridge (UK) and New York 
Radhakrishnan V., Srinivasan G., 1982, Current Science, 
51, 1096 

Romanova M. M., Ustyugova G. V., Koldoba A. V., 

Lovelace R. V. E., 2004, ApJ, 610, 920 
Stairs I. H., Lyne A. G., Shemar S. L., 2000, Nature, 406, 

484 

Strohmayer T. E., Markwardt C. B., Swank J. H., in't Zand 

J., 2003, ApJL, 596, L67 
Ushomirsky G., Cutler C, Bildsten L., 2000, MNRAS, 319, 

902 

Vanderspek R., Morgan E., Crew G., Graziani C, Suzuki 

M., 2005, The Astronomer's Telegram, 516, 1 
Vigelius M., Melatos A., 2008, ArXiv e-prints, 802 
Watts A., Krishnan B., Bildsten L., Schutz B., 2008, ArXiv 
e-prints, 803 

Watts A. L., Patruno A., van der Klis M., 2008, ArXiv 
e-prints, 805 

Watts A. L., Strohmayer T. E., 2006, MNRAS, 373, 769 
Watts A. L., Strohmayer T. E., Markwardt C. B., 2005, 
ApJ, 634, 547 

Wijnands R., 2004, Nuclear Physics B Proceedings Supple- 
ments, 132, 496 

Wijnands R., 2006, in Lowry J. A., ed., Trends in Pulsar 
Research Accretion-Driven Millisecond X-ray Pulsars, pp 
53-+ 

Zimmermann M., Szedenits Jr. E., 1979, Phys. Rev. D, 20, 
351 

This paper has been typeset from a TffX/ DTffX file prepared 
by the author. 



